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Abstract 



A search for second and third generation scalai^ and vector leptoquarks produced in ep 
collisions via the lepton flavour violating processes ep — pX and ep — )■ tX is performed 
by the HI experiment at HERA. The full data sample taken at a centre-of-mass energy 
^/s = 319 GeV is used for the analysis, coiTcsponding to an integrated luminosity of 
245 pb~^ of e~^p and 166 pb~^ of e~p collision data. No evidence for the production of 
such leptoquarks is observed in the HI data. Leptoquarks produced in e^p collisions with 
a coupling strength of A = 0.3 and decaying with the same coupling strength to a muon- 
quark pair or a tau-quark pair are excluded at 95% confidence level up to leptoquark masses 
of 712 GeV and 479 GeV, respectively. 
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1 Introduction 



The ep collisions at HERA provide a unique testing ground to search for new particles coupling 
directly to a lepton and a quark. An example of such particles are leptoquarks (LQs), colour 
triplet bosons which are a generic prediction of grand unified theories [[I], composite models [|2l, 
technicolour [Sj and supersymmetry with i?-parity violation dH. In the Standard Model (SM) 
particle interactions conserve lepton flavour, and if this property is extended to LQ models, any 
such particles produced at HERA would decay exclusively into a quark and a first generation 
lepton, namely an electronQor a neutrino. Dedicated searches have been performed at HERA for 
such leptoquarks, where the SM expectation is dominated by neutral current (NC) and charged 
current (CC) deep inelastic scattering (DIS) background [l5]-|3. 

The introduction of lepton flavour violation (LEV) to leptoquark models would mean the 
processes ep — )■ fiX or ep — rX, mediated by the exchange of a leptoquark, would be ob- 
servable at HERA with final states containing a muon or the decay products of a tau lepton in 
combination with a hadronic system X. Searches for such signatures have been performed at 
HERA and limits on LEV leptoquark production have been derived [I51[8l|9l. In this paper a 
search for LEV phenomena is performed using e^p collision data at a centre-of-mass energy 

= 319 GeV, recorded during the years 1998-2007 by the HI experiment at HERA. The cor- 
responding integrated luminosity of 245 pb~^ for e^p collisions and 166 pb~^ for e~p collisions 
represents an increase in size of the data sample with respect to the previous publication by a 
factor of 3 and 12, respectively. Data collected from 2003 onwards were taken with a longitudi- 
nally polarised lepton beam, with polarisation typically at a level of 35%. The presented results 
supersede those derived in previous searches for lepton flavour violating leptoquarks by the HI 
experiment (HI . 



2 Leptoquark Phenomenology 

The phenomenology of LQs at HERA is discussed in detail elsewhere [21. In the framework of 
the Buchmiiller-Ruckl-Wyler (BRW) effective model [flOl , LQs are classified into 14 types IfTTI 
with respect to the quantum numbers spin J, weak isospin / and chirality C{= L, R). Scalar 
(J = 0) LQs are denoted as Sf and vector (J = 1) LQs are denoted Vf in the follow- 
ing. LQs with identical quantum numbers but different weak hypercharge are distinguished 
using a tilde, for example Vq^ and V^. Some LQs, namely Sq, S^, Vq and Vj^, may de- 
cay to a neutrino-quark pair resulting in the branching fraction for decays into charged leptons 
f3i = Tig/ {Tig + T^^g) =0.5. Slncc neutrino flavours cannot be distinguished with the HI exper- 
iment such final states are not included in this analysis. 

Leptoquarks carry both lepton (L) and baryon (B) quantum numbers, and the fermion 
number F = L + 3B is assumed to be conserved. Leptoquark processes proceed directly via 
s-channel resonant LQ production or indirectly via n-channel virtual LQ exchange. Eor LQ 
masses rriLQ well below y/s, the s-channel production of F = 2 (F = 0) LQs in e~p (e~^p) 

'In this letter the term "electron" is used generically to refer to both electrons and positrons, if not otherwise 
stated. 
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collisions dominates. However, for LQ masses above 319 GeV, both the s and w-channel pro- 
cesses are important such that both e~p and e+p collisions have similar sensitivity to LQs with 
F = 2 and LQs with F = 0. 

The BRW model assumes lepton flavour conservation, although a general extension of this 
model allows for the decay of LQs to final states containing a quark and a lepton of a different 
flavour, that is a muon or tau lepton. Non-zero couplings X^g- to an electron-quark pair and A^^^ 
(Xrqj) to a muon(tau)-quark pair are assumed. The indices i and j represent quark generation 
indices, such that X^q^ denotes the coupling of an electron to a quark of generation i, and X^g^ is 
the coupling of the outgoing lepton (where £ = /i or r) to a quark of generation j. An overview 
of this extended model for the LQ coupling to u and d quarks is provided elsewhere [H]. 

Events with LQs are generated using the LEGO [[T2l event generator with the CTEQ5L 
parametrisation [[T3l of the parton distribution functions of the proton. The LQ signal expecta- 
tion is calculated as a function of the LQ type, mass, coupling constant and the branching ratio 
/3 to a given charged lepton flavour, where: 

{3 = (3,x P,^y with P,,y = -^^^^^ and T,, = m^^Xl x '""^"'^^ 

i /.(r)g + i eg [ 24^ VCCtOr LQ 

where T^g denotes the partial LQ decay width for the decay to a lepton i = e, ii,t and a quark 
q. In order to avoid the need to generate many Monte Carlo (MC) samples at each leptoquark 
mass, coupling and branching ratio, a weighting technique is used to provide predictions across 
the full range of LQ production parameters [l8] . 

Leptoquarks with couplings to the first and the second lepton generation may decay to a 
muon and a quark, leading to an event topology with an isolated high transverse momentum Pt 
muon back-to-back to a hadronic system in the transverse plane. Leptoquarks with couplings 
to the first and the third lepton generation may decay to a tau and a quark. Tau leptons are 
identified in this analysis using the muonic and one-prong hadronic decays of the tau. In both 
cases, the tau decay results in missing transverse momentum in the event due to the escaping 
neutrinos. Previous LEV leptoquark analyses also examined r — > eX decays and three-prong 
hadronic tau decays [[8l, for which the background from SM processes is large [[T4l . Given the 
increase in data luminosity with respect to the previous publication, a correspondingly large 
increase in the SM background is observed, which limits the sensitivity of these decay channels 
and they are therefore not included in the presented analysis. 



3 Standard Model Background Processes 

Several SM processes may mimic the LQ signal. The main SM background contribution is 
from photoproduction events, in which a hadron is wrongly identified as a muon or a narrow 
hadronic jet fakes the signature of the hadronic tau decay. Similarly, the scattered electron in 
NC DIS events may also be misinterpreted as the one-prong hadronic tau decay jet. Smaller 
SM background contributions arise from events exhibiting intrinsic missing transverse momen- 
tum (for example CC DIS), events containing high Pt leptons (such as lepton pair production. 
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particularly inelastic muon-pair events if one muon is unidentified) or events with both of these 
features (real W production with leptonic decay). 

The RAPGAP ifTSl event generator, which implements the Bom level, QCD Compton and 
boson-gluon fusion matrix elements, is used to model inclusive NC DIS events. The QED ra- 
diative effects arising from real photon emission from both the incoming and outgoing electrons 
are simulated using the HERACLES [161 program. Direct and resolved photoproduction of jets 
and prompt photon production are simulated using the PYTHIA [17] event generator, which is 
based on Bom level scattering matrix elements. In RAPGAP and PYTHIA, jet production from 
higher order QCD radiation is simulated using leading logarithmic parton showers and hadroni- 
sation is modelled with Lund string fragmentation [[T8l . Inclusive CC DIS events are simulated 
using the DJANGOH |[T9ll program, which includes first order leptonic QED radiative correc- 
tions based on HERACLES. The production of two or more jets in DJANGOH is accounted 
for using the colour dipole model ll20l . The leading order MC prediction of processes with two 
or more high transverse momentum jets in NC DIS, CC DIS and photoproduction is scaled by 
a factor of 1.2 to account for the incomplete description of higher orders in the MC genera- 
tors ll2Tll22]| . Contributions arising from the production of single W bosons and multi-lepton 
events are modelled using the EPVEC ll23l and GRAPE [|24ll event generators, respectively. The 
uncertainties on the SM background predictions are described in section [6l 

Generated events are passed through a GEANT [|25l based simulation of the HI appara- 
tus, which takes into account the mnning conditions of the data taking. Simulated events are 
reconstmcted and analysed using the same program chain as is used for the data. 

4 Experimental Conditions 

A detailed description of the HI experiment can be found elsewhere [|26l . Only the detec- 
tor components relevant to this analysis are briefly described here. A right-handed Cartesian 
coordinate system is used with the origin at the nominal primary ep interaction vertex. The 
proton beam direction defines the positive z axis (forward direction). The polar angle 9 and 
the transverse momenta Pt of all particles are defined with respect to this axis. The azimuthal 
angle (p defines the particle direction in the transverse plane. The pseudorapidity is defined as 
?7 = — In tan |. 

The Liquid Argon (LAr) calorimeter ['271 covers the polar angle range 4° < 6' < 154° with 
full azimuthal acceptance. The energies of electromagnetic showers are measured in the LAr 
with a precision of a{E)/E ^ 11%/v/^/GeV© 1% and hadronic energy depositions with 
a{E)/E ~ 50%/^/E/GeM © 2%, as determined in test beam measurements [|28l|29l. A lead- 
scintillating fibre calorimeter (SpaCal) [l30l covering the backward region 153° < 9 < 178° 
completes the measurement of charged and neutral particles. For electrons a relative energy 
resolution of a{E)/E ~ 7%/ ^/E/ GeV © 1% is reached, as determined in test beam measure- 
ments [13I]|- The central (20° < 9 < 160°) and forward (7° < 6* < 25°) inner tracking detectors 
are used to measure charged particle trajectories and to reconstmct the interaction vertex. The 
measured trajectories fitted to the interaction vertex are referred to as tracks in the following. 
The LAr calorimeter and inner tracking detectors are enclosed in a superconducting magnetic 
coil with a field strength of 1.16 T. From the curvature of charged particle trajectories in the 
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magnetic field, tlie central tracking system provides transverse momentum measurements with 
a resolution of ap^/PT = 0.005PT/GeV © 0.015 Il32l . The return yoke of the magnetic coil is 
the outermost part of the detector and is equipped with streamer tubes forming the central muon 
detector (4° < 6* < 171°). In the very forward region of the detector (3° < 9 < 17°) a set of 
drift chambers detects muons and measures their momenta using an iron toroidal magnet. The 
luminosity is determined from the rate of the Bethe-Heitler process ep — )■ ep7, measured using 
a photon detector located close to the beam pipe at 2; = —103 m, in the backward direction. 

Lepton flavour violating processes usually exhibit an imbalance in the measured calorimet- 
ric transverse momentum, P^^^", due to either the presence of a minimally ionising muon in 
HX final states or the escaping neutrino(s) from tau decays in rX events. The LAr calorime- 
ter provides the main trigger in this analysis. The trigger efficiency is about 60% for events 
with a transverse momentum imbalance measured in the calorimeter of 12 GeV, rising to about 
90% for an imbalance of 25 GeV [|33l . Events are also triggered by hadronic jets in the LAr 
calorimeter, with a trigger efficiency above 95% for a jet transverse momentum P^^ > 20 GeV 
and almost 100% for P^^ > 25 GeV Il34l . For di-jet events with a scalar sum of the transverse 
energy in the event Et > 30 GeV, the trigger efficiency is greater than 98% Il35l . 

In order to remove events induced by cosmic rays and other non-ep background, the event 
vertex is required to be reconstructed within ±35 cm in z of the average nominal interaction 
point. In addition, topological filters and timing vetoes are applied. 



5 Particle Identification and Event Selection 

Electromagnetic particle (electron and photon) candidates are identified as compact and iso- 
lated clusters of energy in the electromagnetic part of the LAr calorimeter. Electron candidates 
are defined as electromagnetic particle candidates with an associated track. Identification of 
muon candidates is based on a track in the inner tracking detectors, associated to a signal in 
the muon system. Tracks and calorimeter deposits not identified as originating from isolated 
electromagnetic particles or muons are combined into cluster-track objects to reconstruct the 
hadronic final state [|36l . Jets are reconstructed from these objects using an inclusive kx algo- 
rithm Il37ll38l with a minimum P^ of 4 GeV and a distance parameter R = 1.0. The missing 
transverse momentum P^^^^, which may indicate the presence of neutrinos in the final state, 
is derived from all reconstructed particles in the event. The LQ kinematics are reconstructed 
using the double angle method ll39l . The direction of the detected lepton and the hadronic fi- 
nal state are used to reconstruct the Bjorken scaling variable x and subsequently the LQ mass 

5.1 Search for second generation leptoquarks 

An initial sample of events with muons and jets is selected by requiring at least one P^ > 8 GeV 
muon in the polar angular range 10° < 9^ < 120° and at least one jet. In addition, P^^'° is 
required to be greater than 12 GeV. After this selection, 996 events are observed in the data. 
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in good agreement with the SM prediction of 978 ± 187, where the uncertainty includes the 
statistical and systematic errors (see section |6]). 

Events with at least one isolated muon are then selected, which is done by requiring the 
angular distance, D = ^ (Ar/)^ + (A0)2, of the muon to the nearest track and to the nearest 
jet to be greater than 0.5 and 1.0, respectively. In addition, an isolated muon may have no 
more than 5 GeV deposited in the LAr calorimeter within a cylinder centred on the muon track 
direction of radius 35 cm (75 cm) in the electromagnetic (hadronic) section. The muon isolation 
requirements reduce the number of selected data events to 220, compared to a SM prediction of 
218 ±48. 

The NC DIS background is further suppressed by increasing the cut on the calorimetric 
momentum imbalance to P^^^° > 25 GeV, which implicitly increases the minimum muon trans- 
verse momentum, and by rejecting events with identified isolated electrons. To reduce the 
muon-pair SM background, exactly one isolated muon is required, as expected in LFV LQ sig- 
nal events. The back-to-back event topology in the azimuthal plane is also exploited to remove 
the SM background and the difference between the azimuthal angle of the hadronic system and 
the muon A0^_x is required to be greater than 170°. As the majority of the energy deposited 
in the calorimeter is due to the hadronic final state, signal events tend to exhibit an azimuthal 
imbalance when considering the calorimeter measurement alone. Therefore, a requirement of 
Vap/Vp < 0.3 is also employed, where Kp/Vp is defined as the ratio of the anti-parallel to 
parallel projections of all energy deposits in the calorimeter with respect to the direction of 
P^^^° ll40l . After these selection cuts, the data sample is reduced to 6 events, compared to a SM 
prediction of 7.5 ± 1.8. 

To exploit the longitudinal balance of the event, a requirement on the sum of the energy 
and longitudinal momentum of all detected particles i in the event Tji{E^ — PI) > 40 GeV is 
applied. In the case of signal events this quantity is expected to be around 2E^ = 55.2 GeV, 
where -Eg is the electron beam energy. However, for the remaining SM background after the 
above event selection the scattered electron or some other backward going final state particle is 
typically undetected, resulting in significantly lower values of — P^). In order to improve 
the resolution, which is poor for very high Pt muons due to the small curvature of the track, the 
transverse momentum of the muon is calculated from the hadronic system, = —P^ and the 
muon track direction is used to reconstruct the longitudinal component P^ and energy E^^ used 
in the (E - P,) sum iHTIl . 

The cut on T.i{E^ — PI) removes five of the remaining data events, so that one event is 
observed in the final selection of the analysis of jiX final states, which compares well to the 
SM prediction of 2.0 ± 0.4, where the largest contribution comes from muon-pair events. The 
presented analysis has a lower background contamination and an improved selection efficiency 
with respect to the previous HI publication. The selection efficiency typically ranges from 75% 
for LQs masses of around 150 GeV to 65% for LQ masses above 300 GeV, representing an 
improvement of an additional 15 — 25% with respect to the previous publication [|8l. 

5.2 Search for third generation leptoquarks 

In the search for third generation leptoquarks, tau leptons are identified in this analysis using 
the muonic and one-prong hadronic decays of the tau. 
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Muonic tau decays r — t- jiu^Vj. result in final states similar to the high Pt muon signatures 
described in section 15.11 The same selection cuts are therefore applied in this channel. To 
account for possible effects due to different muon kinematics resulting from the tau decay, 
the selection efficiency was studied in an LFV MC signal sample with a tX final state and a 
subsequent muonic tau decay. The selection efficiency in this channel is up to 70% at leptoquark 
masses of around 150 GeV and about 55% for LQ masses above 300 GeV, which represents a 
similar level of improvement with respect to the previous publication ^ as observed in the 
second generation search described in section 15.11 

The one-prong hadronic decay of the tau leads to a high Pt, narrow "pencil-like" jet, so that 
the typical LFV signal event topology is a di-jet event. An initial event sample for the analysis 
of this decay channel is formed by selecting events with at least two jets in the polar angle 
range 5° < ^j"* < 175° and with Pj?*^ > 20 GeV and P^^"^ > 15 GeV. This results in a large 
di-jet sample of approximately 2.2 ■ 10^ events, which is consistent with the SM prediciton of 
(2.6 ± 0.5) ■ 10^, where the main contribution comes from photoproduction. 

A tau jet is characterised by a narrow energy deposit in the calorimeter and a low track 
multiplicity within the identification cone of the jet. Tau jet candidates are identified in the di- 
jet sample, where the candidate is required to be in the polar angle range 20° < < 120° and 
has a maximum jet radius i?jct of 0.12 [42]. The jet radius is used as a measure of the coUimation 
of the jet and is calculated as: i?jct = -E'h\/^'7(jet, h)'^ + A0(jet, hy, where i^jct is the 

total jet energy and the sum runs over all jet daughter hadronic final state particles of energy 
Eh- At least one track with P-r larger than 2 GeV not associated with an identified electron or 
muon is required within the jet radius of the tau jet candidate. Approximately 3 ■ lO'' tau jet 
candidates are identified in the di-jet sample. 

The undetected neutrinos from tau lepton decays result in an overall Pt imbalance and 
therefore a minimum missing transverse momentum P^^^^ > 12 GeV is required. Events with 
only one tau jet are then selected, which is required to be isolated from tracks and other jets 
in the event by a distance D > 1.0. A track multiplicity of one is required in a cone of radius 
R = 1.0 around the jet axis. Tau jets with additional track segments not fitted to the event ver- 
tex within a cone of radius R = 0.3 around the jet axis are also rejected [43]. To reject purely 
electromagnetic jets, a maximum of 90% of the jet energy may be recorded in the electromag- 
netic part of the calorimeter. The resulting selection contains 104 data events compared to a SM 
prediction of 116 ± 16. 

Further cuts are then applied to reduce the remaining SM background. The hadronic trans- 
verse momentum P^ is required to be larger than 30 GeV and the acoplanarity between the tau 
jet and X system in the transverse plane Acpr^x is required to be greater than 160°. Note that 
for the analysis of the hadronic tau decay channel, the tau jet is subtracted from the inclusive 
hadronic final state to obtain the four- vector of the remaining hadronic system X. Analogous to 
the muon channel, a cut of T,i(E'^ — PI) > 40 GeV is applied to exploit the longitudinal balance 
of the event. Similarly to the muon channel, only the direction of the tau jet is used in the sum, 
and the transverse momentum of the X system is employed to determine the tau jet four- vector. 
Electrons entering inactive regions of the electromagnetic section of the LAr calorimeter may 
fake the tau jet signature and therefore these regions are excluded from the analysis [|4TI . 

In the analysis of rX final states where the tau lepton decays hadronically, 6 events are 
observed in the data, in good agreement with the SM prediction of 8.2 ± 1.1, where the main 



9 



SM contribution is from remaining NC DIS events. The selection efficiency in the one-prong 
hadronic tau decay channel ranges from 18% for leptoquark masses in the range 150-200 GeV 
to 12% for masses above 300 GeV. 



6 Systematic Uncertainties 

The following experimental systematic uncertainties are considered in the search for second 
generation leptoquarks: the scale uncertainty on the transverse momentum of high Pt muons 
is 2.5% and the uncertainty on the muon polar angle measurement 3 mrad p4]; the muon 
identification efficiency has an error of 5% in the region 9^^ > 12.5° and 15% in the forward 
region [|43l ; the hadronic energy scale is known within 2% and the uncertainty on the hadronic 
polar angle measurement is 10 mrad [|34l . In the search for third generation leptoquarks an 
uncertainty on the description of the jet radius _Rjct is included in the analysis by varying the 
cut value of 0.12 by 10%. All other experimental systematic uncertainties in the tau channel are 
included in the model uncertainties described below [43]. In both searches, the uncertainty on 
the trigger efficiency is 2-3% and the uncertainty on the luminosity measurement is 3%. 

The effects of these systematic uncertainties on the signal and the expected SM background 
are evaluated by shifting the relevant quantities in the MC simulation by their uncertainty and 
adding all resulting variations in quadrature. 

Additional model uncertainties are attributed to the normalisation uncertainties in the analy- 
sis phase space of the SM MC generators described in section |3l These model uncertainties are 
estimated from control analyses in an extended phase space relevant to the search signature |P^3l . 
In the analysis of fiX final states, the contributions from RAPGAP (NC DIS), PYTHIA (photo- 
production) and GRAPE (lepton-pair production) are each attributed a systematic error of 30%, 
which is increased to 50% for the period 1998-2000 [Al}. The contribution from DJANGOH 
(CC DIS) in events with isolated muons is attributed an uncertainty of 50%. In the analysis 
of tX final states where the tau lepton decays hadronically, the contribution from RAPGAP, 
PYTHIA, DJANGOH and GRAPE are attributed systematic uncertainties of 15%, 20%, 20% 
and 30%, respectively. The theoretical uncertainty of 15% is used for all predicted contributions 
from EPVEC (W production) [l23l . 

The total error on the SM prediction is determined by adding the MC statistical error to the 
effects of all model and experimental systematic uncertainties in quadrature. 

The main theoretical uncertainty on the signal cross section originates from the parton den- 
sities. This uncertainty is estimated to be 5% for LQs coupling to up-type quarks and varies 
between 7% at low masses and 30% at masses around 290 GeV for LQs coupling to down-type 
quarks 0. 

7 Results 

The observed number of events is in agreement with the SM prediction and therefore no evi- 
dence for LEV is found. The reconstructed leptoquark mass in the search for ep — >■ fxX and 
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ep — )■ tX events is shown in figured] compared to the SM prediction and an example LQ signal 
with arbitrary normalisation. 

In the absence of a signal, the results of the search are interpreted in terms of exclusion limits 
on the mass and the coupling of LQs that may mediate LFV. The LQ production mechanism 
at HERA involves non-zero coupling to the first generation fermions Agg > 0. For the LFV 
leptoquark decay, it is assumed that only one of the couplings A^g and A^^ is non-zero and that 
Aeg = A^g(AT-g), which results in I^lfv = 0-5. A modified frequentist method with a likelihood 
ratio as the test statistic is used to combine the individual data sets and the ep — > tX search 
channels P31 . The lepton beam polarisation enters the limit calculation for the 2003-2007 data. 

Figures |2]and|3]display the 95% confidence level (CL) upper limits on the coupling X^q and 
\rq of all 14 LQ types to a muon-quark pair and a tau-quark pair, respectively, as a function of 
the mass of the LQ leading to LFV in ep collisions. Only first generation quarks are considered 
in these limits. The limits are most stringent at low LQ masses with values of C(10~'^) at 
= 100 GeV The limits corresponding to LQs coupling to a n quark are more stringent 
than those corresponding to LQs coupling to the d quark only, as expected from the larger 
u quark density in the proton. Corresponding to the steeply falling parton density function 
for high values of x, the LQ production cross section decreases rapidly and exclusion limits 
are less stringent towards higher LQ masses. For LQ mass values near the kinematic limit of 
319 GeV, the limit corresponding to a resonantly produced LQ turns smoothly into a limit on 
the virtual effects of both an off-shell s-channel LQ process and a n-channel LQ exchange. At 
masses ttt-lq ^ \fs the two processes contract to an effective four-fermion interaction, where 
the cross section is proportional to (AegA^(^)g/mLQ)^. For a coupling A of electromagnetic 
strength, where A = ^/4^ra^ = 0.3, LFV leptoquarks produced in ep collisions decaying to a 
muon-quark or a tau-quark pair are excluded at 95% confidence level up to leptoquark masses 
of 712 GeV and 479 GeV, respectively. 

The limits on \eq = Xfi{T)q in the region ^tilq ^ are transformed into a limit on the 
value Aeg- A^(T-)g^ /ttilq and shown in tables [T] and [2] for F = LQs and in tables [3] and |4] for 
F = 2 LQs. For each LQ type, the limit is calculated for the hypothesis of a process with only 
the quarks of flavours i and j involved. With respect to quark flavours, the selection criteria 
described in sections 15.11 and 15.21 are inclusive since no flavour tagging of the hadronic jet is 
used. Nevertheless, the sensitivity of the analysis to different quark flavours varies due to the 
parton content of the proton and the presence of the n-channel exchange. Leptoquark couplings 
to the top quark are not considered in these limits. 

The HI limits may be compared with constraints from low energy experiments, based on the 
non-observation of LFV in muon scattering and rare decays of mesons and leptons [46.1 . The 
interpretation of these results in terms of leptoquark exchange and limits on ^eq^ii{T)qj I'^i^q BTl 
are also shown in the tables . Superior limits are observed by H 1 in the search for third generation 
leptoquarks, compared to limits from B ^ re decays, as well as in a few unique channels. 

At hadron colliders, LQs are mainly produced in pairs independently of A, and therefore 
searches cannot constrain the LFV couplings. Lower mass limits by the CMS experiment on 
second generation scalar leptoquarks extend up to 394 GeV [48] for a branching ratio (3 = 1. 
Third generation scalar (vector) leptoquarks are ruled out below 247 GeV by the D0 experi- 
ment [HI (317 GeV by the CDF experiment [EOl) for /3 = 1. For /3 = 0.5, which is a more 
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appropriate value for a comparison to the production of LQs at HERA, the above second gen- 
eration search rules out leptoquark masses below around 300 GeV, at which mass this analysis 
rules out such scalar LQs with couplings in the range A = 0.2-0.3. 



8 Conclusion 

A search for lepton flavour violation processes induced by leptoquarks in ep collisions at a 
centre-of-mass energy of 319 GeV with the HI experiment at HERA is presented. No signal 
for the LEV processes ep — )■ jjiX or ep — > tX is observed and assuming a coupling strength of 
A = 0.3, leptoquarks mediating lepton flavour violation are ruled out up to masses of 712 GeV 
and 479 GeV, respectively. The new HI limits extend beyond the domain in LQ mass excluded 
by previous searches at HERA. Additionally, the HI limits remain competitive in certain chan- 
nels with those from low energy experiments and for large values of the couplings exclude 
leptoquark masses beyond the current limits from hadron colliders. 
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ep^ fiX HI F = 


Upper exclusion limits on \pq^\^^q./m\(^ (TeV^) 
for lepton flavour violating leptoquarks at 95% CL 
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Table 1: Limits at 95% CL on Xeq-X^g./mlQ for F = leptoquarks (bold). The fermion pairs 
considered in the analysis coupling to each LQ type are indicated in the column headings. The 
5*^2 ^1 LQs couple to both w-type {U) and d-type {D) quarks [|TOl . The cases marked 
with refer to scenarios involving a top quark. Combinations of i and j shown in the first 
column denote the quark generation coupling to the electron and muon respectively. In each 
cell the first two rows show the process providing the most stringent limit from low energy 
experiments. Highlighted HI limits are more stringent than those from the corresponding low 
energy experiment. 
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ep^rX HI F = 



Upper exclusion limits on Xeq,^XrqJm1q (TeV ^) 
for lepton flavour violating leptoquarks at 95% CL 
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Table 2: Limits at 95% CL on Xeg-XTg^ / fnl^ for F = leptoquarks (bold). The fermion pairs 
considered in the analysis coupling to each LQ type are indicated in the column headings. The 
5*^2 ^1 LQs couple to both w-type {U) and rf-type {D) quarks [|TOl . The cases marked 
with refer to scenarios involving a top quark. Combinations of i and j shown in the first 
column denote the quark generation coupling to the electron and tau lepton respectively. In 
each cell the first two rows show the process providing the most stringent limit from low energy 
experiments. Highlighted HI limits are more stringent than those from the corresponding low 
energy experiment. 
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ep^ fiX HI F = 2 


Upper exclusion limits on \pq^\^^q./m\(^ (TeV^) 
for lepton flavour violating leptoquarks at 95% CL 
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Table 3: Limits at 95% CL on Xeq-X^g./mlQ for F = 2 leptoquarks (bold). The fermion pairs 
considered in the analysis coupling to each LQ type are indicated in the column headings. The 
and Vf^2 LQs couple to both w-type (U) and d-type (D) quarks [|TOl . The cases marked 
with refer to scenarios involving a top quark. Combinations of i and j shown in the first 
column denote the quark generation coupling to the electron and muon respectively. In each 
cell the first two rows show the process providing the most stringent limit from low energy 
experiments. Highlighted HI limits are comparable to those from the corresponding low energy 
experiment. 
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Table 4: Limits at 95% CL on Xeg-XTg^ / fnl^ for F = 2 leptoquarks (bold). The fermion pairs 
considered in the analysis coupling to each LQ type are indicated in the column headings. The 
and Vf^2 LQs couple to both w-type (U) and d-type (D) quarks [|TOl . The cases marked 
with refer to scenarios involving a top quark. Combinations of i and j shown in the first 
column denote the quark generation coupling to the electron and tau lepton respectively. In 
each cell the first two rows show the process providing the most stringent limit from low energy 
experiments. Highlighted HI limits are more stringent than those from the corresponding low 
energy experiment. 
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Figure 1: The reconstructed leptoquark mass in the search for (a) ep jiX and (b) ep rX 
events. The data are the points and the total uncertainty on the SM expectation (open histogram) 
is given by the shaded band. The dashed histogram indicates the LQ signal with arbitrary 
normalisation for a leptoquark mass of 150 GeV. 
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Figure 2: Exclusion limits on the coupling constants A^g = Xeq as a function of the leptoquark 
mass Mlq for (a) scalar LQs with F = 0, (b) vector LQs with F = 0, (c) scalar LQs with 
F = 2 and (d) vector LQs with F = 2. Regions above the lines are excluded at 95% CL. 
The notation qi indicates that only processes involving first generation quarks are considered. 
The parentheses after the LQ name indicate the fermion pairs coupling to the LQ, where pairs 
involving anti-quarks are not shown. 
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Figure 3: Exclusion limits on the coupling constants \rq = Xeq as a function of the leptoquark 
mass Mlq for (a) scalar LQs with F = 0, (b) vector LQs with F = 0, (c) scalar LQs with 
F = 2 and (d) vector LQs with F = 2. Regions above the lines are excluded at 95% CL. The 
notation qi indicates that only processes involving first generation quarks are considered. The 
parentheses after the LQ name indicate the fermion pairs coupling to the LQ; pairs involving 
anti-quarks are not shown. 
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